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Abstract 

It is shown that the wounded nucleon model describes very well 
the recent PHOBOS data on particle production in D-Au collisions at 
200 GeV. Contribution to particle production from a single wounded 
nucleon is determined. A two-component model is formulated and 
shown to account for most of the important features of the data. 

1 Introduction 

The model of wounded nucleons, proposed almost 30 years ago pQ, shows a 
remarkable survival capacity: it is still being used in analysis of data 0] 
and the very concept of a "wounded" nucleon (called now a "participant" 1 ) 
became one of the basic tools in description and interpretation of the heavy 
ion experiments. 

In its original form, the model proposes that the particle production in 
a nucleus-nucleus collision can be represented as a superposition of indepen- 
dent contributions from the wounded nucleons in the projectile and in the 

*Address: Reymonta 4, 30-059 Krakow, Poland; e-mail:bialas@th. if. uj.edu. pi; 

1 According to the definition given in £Q, the wounded nucleon is the one which un- 
derwent at least one inelastic collision. We stick to this name because we think that the 
name "participant" is incorrect: It should rather refer to nucleons which underwent any 
(either elastic or inelastic) collision. 
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target. Consequently, the density of particles in a collision of nuclei of nuclear 
numbers A and B is given by 



N j AB = w A F A (y) + w B F B (y) = 
dy 

= \{w A + w B )[F A {y) + F B {y)\ + ±(w A - w B )[F A {y) - F B {y)]lwl (1) 

where w A and w B are the numbers of the wounded nucleons in A and B, y 
is the rapidity in the cm. system of the collision and F A (y) is a contribution 
from a single wounded nucleon in A. Similarly, F B (y) is the contribution 
from a single wounded nucleon in B. The model requires 

F B (y) = F A {-y)lwwl (2) 

but it will be convenient to keep the more general formalism. 

Recently, the pseudorapidity distribution dN/dr] of particles produced in 
d-Au collisions was measured by PHOBOS and BRAHMS collaborations at 
RHIC in a wide range of available phase-space and for various centralities 
jHEJ. In the present paper we use the wounded nucleon model to analyze 
the data reported by PHOBOS 0|. 

We find that the model gives a good description of the data, with the 
condition (??) being well satisfied, except at rapidities close to the maximal 
values. This observation allows us to determine from the data the contribu- 
tion F(r]) from a single wounded nucleon. 

Two novel features emerge from this analysis. It turns out that (i) F(rf) 
is not confined to the hemisphere corresponding to the wounded nucleon 
in question but rather extends over all available rapidity (except possibly 
close to the boundary); Moreover, (ii) F{rf) shows a distinct two-component 
structure. 

We have shown that these observations can be understood in a recently 
proposed mechanism 0, describing the particle production as a two-step 
process: (i) multiple colour exchanges between partons from projectile and 
target, followed by (ii) particle emission from colour sources created in the 
first step. 

In the next section we show that the wounded nucleon model describes 
correctly the data from PHOBOS 4J. Determination of the contributions 
F Au (rj) and Fnirj) from the wounded nucleons is presented in Section 3 where 
also their structure is discussed. In Section 4 we discuss a possible explana- 
tion of these findings. Our conclusions are listed in the last section. 
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2 Wounded nucleons in Deuteron-Gold colli- 



sions 

The most direct way to test the relation (??) is to construct the symmetric 
and antisymmetric components of the particle density: 



drj 



drj 



(3) 



In Figures 1 and 2 these two quantities are plotted versus pseudorapidity 
for various centralities measured in |4. To compare with the model, we 
constructed the averages 



E>£±4 c) ]/2 



-I2w 



(4) 



The model predicts [c.f. (??)] 



^ - < > I3w 



(5) 



The R.H.S of (??) is shown in Figs 1 and 2 as shaded areas (expressing the 
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Figure 1: Antisymmetric part of the deuteron D-Au inclusive cross-section 
compared with predictions of the wounded nucleon model. 



inacuracies in determination of wa u and wd). One sees that the agreement 
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Figure 2: Symmetric part of the deuteron D-Au inclusive cross-section com- 
pared with predictions of the wounded nucleon model. 



is rather satisfactory (except in the regions close to the maximal allowed 
rapidity) . 

We thus conclude that the model describes correctly the available data. 

Using data on nucleon-nucleon collisions another, more demanding, test 
of the model is possible. 

Indeed, one sees immediately from (??) that for the nucleon-nucleon col- 
lision we have 

= F N (y) + F N (-y)lw2 (6) 

and thus for the ratio 

dN AB /dy 

RAB{y) = ~7T7 — jt 1w3 7 
dN NN dy 



one obtains 



Rab{ V ) = \{w A + w B ) + \{w A - w B) ^J ) + p^ lwA (8) 
The first immediate consequence is 

Rab(v = 0) = ~(w A + w B )lw5 (9) 



4 



implying that the value of the ratio Rab at mid-rapidity is fully determined 
by the number of wounded nucleons and entirely independent of the shape 
of the function F(y). 

Fig. 3 shows Rd-Au(0) plotted versus (wa u + wn)/2, as measured by the 
Phobos collaboration [3] 2 . One sees that the data are indeed in excellent 
agreement with (??). 




Figure 3: Particle production in the central region compared with the pre- 
dictions of the wounded nucleon model. 

The next step is to verify if the model gives an adequate description of 
data at y ^ 0. To this end we propose to study the quantity 

D D -Au{v) = dNl l~ Au - R d ^au(v = 0) dN f N lww (10) 
aq arj 

which, according to (??) and (??) should obey 

D Au -D{y) = ^(w Au - w D )<5>{y)lw8 (11) 

where 



^(y)=F A u(y)-F D (y).lw9 (12) 

In Fig. 4 Dd_au(v) i s plotted for various centralities, as measured in the 
PHOBOS experiment. 



2 The numerator of -Rd-^^O) was taken from the numerical data given in [1]. The 
denominator was read off from the Figure lb of the same paper. 
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Figure 4: Comparison of the relation (10) with the predictions of the wounded 
nucleon model. 



To verify (??) we again construct the "average over centralities'' 



£>ii-™g } ]/2 



< m >= „r ( / ( C ; ,> 10 ( l3 ) 



where c denotes the centrality, as measured by PHOBOS. The product 
\{wau — wd) < > 1S shown in Fig. 4 as shaded areas. One sees 

good agreement with the measured values of Da u -d(v) m the deuteron hemi- 
sphere 3 . There are deviations in the Au hemisphere for the most central col- 
lisions. They may be either genuine -hitherto unexplained- deviations from 
the model, or simply represent an additional contribution to particle produc- 
tion from the secondary interactions inside the nucleus. More work is needed 
to clarify this feature. 

3 For maximal centrality, the approximate linear dependence on r\ was observed for 
Rd-Au(u) already in We thanks W.Busza for calling our attention to this observations 
which triggered our interest in the subject. 
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3 Particle emission from a single wounded 
nucleon 

From (??) and (??) we deduce that the contribution from a single wounded 
nucleon can be expressed as 

F(±T]) = -[< $+(77) > ± < $-(77) >]lr2 (14) 

The functions F(±rj) and < $ ± (?7) >) are shown in Fig. 5. 
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Figure 5: Particle production from a single wounded nucleon. Symmetrized 
particle densities. 

Similarly, using (??) and (??) and taking into account that the symme- 
try relation (??) is well satisfied by the data, one can express = F D (r]) ps 
Fau(-v), m terms of cr pp (r]) = F D (r}) + F Au (<q) and $(77) = F Au (<q) - F D {rj): 

Fd^^pp + Hv)}] F Au = ^[a pp -$(r])}lrl (15) 

In Fig. 6 Fc{rj) and F Au (r)) are shown in together with (Jp P {rj) and < $(r/) >. 

One sees that, except for tiny details, the results of both figures are 
very similar to each other (if one excludes the regions close to the maximal 
rapidities). We thus conclude that both methods of analysis lead to the same 
picture. 

Three striking features are to be noted: 
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Figure 6: Particle production from a single wounded nucleon. Unsym- 
metrized particle densities. 

(a) One sees that -contrary to naive expectations- the contribution from a 
wounded nucleon extends far beyond its own hemishere, covering practically 
the full rapidity interval except about 1.5 units from both ends (where the 
energy conservation effects and the intranuclear cascade are expected to give 
corrections to the model in any case) 4 . 

(b) Another -hitherto unexpected- observation is the very simple linear 
dependence on 77 of < $(r/) > and of < <$>~(rj) > which can be well approxi- 
mated by a straight line with the slope of about 1/3. 

(c) There is a dramatic difference between the rapidity dependence of the 
symmetric and antisymmetric part of F(i]). 

In the next section will shall discuss the consequences of these observa- 
tions for the mechanism of particle production. 

4 A possible interpretation 

The observations made in the previous section allowed us to determine the 
contribution of one wounded nucleon to particle production. We have thus 
obtained a qualitatively new information on mechanism of the inelastic nucleon- 
nucleon collisions. A possible interpretation of this result is presented below. 

One sees explicitely that the model does not work in the Au fragmentation region, 
where Fd{t]) turns out negative. 
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The striking difference between the measured symmetric and antisym- 
metric part of F(y), seen in Figs 5 and 6, suggests that F(y) (which is the 
sum of its symmetric and antisymmetric parts) may consists of two compo- 
nents of different origin. A natural possibility is to identify one component 
with particle emission from the valence part of the nucleon and another one 
with emission from the gluon cloud. 

We thus write 



where we qualitatively expect the gluonic contribution to dominate the sym- 
metric part of the spectrum, while its asymmetric part is generated by the 
valence contribution. 

To discuss this concept in more detail, we shall use the model proposed 
recently in |Hj , thus accepting that particle production proceeds in two-steps 
(i) the multi-gluon colour exchanges between partons of the projectile and 
of the target and (ii) the following emission of particle clusters from colour 
sources or strings. 

To illustrate the consequences of this idea we shall assume that particle 
emission from a colour source follows the general features of the bremsstrahlung 
mechanism [TO], or -equivalently [H]- the string model 0|H]. Consider a colour 
source moving to the right. It will emit clusters 5 , approximately uniformly 
in rapidity, until it is neutralized by one of the partons of the target. Thus 
the density of the emitted clusters, p(y;y + ,y~), is confined to the rapidity 
region between the rapidity of the source (y + ) and the rapidity (y~) of this 
parton from the target which neutralized the source 6 . 

Consequenly, the observed distribution of clusters is 



where H^- V,9 \y + ) represent the distributions of the emitting sources (valence 
and gluons) and h(y~) the distribution of the partons in the target normalized 
to 1. Assuming that the parton distribution is dominated by gluons we 

5 It is well known that most of the observed particles are decay products of resonances 
or "clusters" [TT]. 

6 In the string language these are rapidities of the two ends of the string. 



dN _ dNW dN^) 
dy dy dy 



lil 



(16) 




(17) 
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furthermore obtain 



%-) = H^(-y-)/ [ Y dyH<*\y)l2i (18) 

J-Y 

This formula immediately implies that the contribution from gluon sources 
is symmetric with respect to y (provided p(y;y + ,y~) is symmetric, as ex- 
pected). Consequently, the contribution to the antisymmetric part of the 
distribution comes solely from the valence sources. 

To illustrate other consequences of (??), we shall first consider a radically 
simplified picture, taking p(y;y + ,y~) = p for y~ < y < y + , and H^ 9 \y) = 
between — Y and Y, where p and H^ 9 ' are constants. The result is 

Since the distribution H^(y + ) of the valence part is confined to the region 
close to maximal rapidity, say y+ > Y*, the integral Jy dy + H^ v \y + ) equals 
1, for y < Y*. Consequently, for y < Y* we have 

i.e. the linear dependence on y. 

This simple exercise shows that the (observed in data) linear dependence 
of the antisymmetric part of the distribution is a direct consequence of the 
flat distribution of gluons and of the emitted clusters. 

The symmetric part is dominated by contribution from gluon sources 
which increases linearly with the total available rapidity. The quadratic 
dependence on y is also -at least qualitatively- not inconsistent with the data 
jS]. Thus we feel that we may be indeed on the right track. 

It is clear that important refinements to this simple example are necessary 
to obtain a more precise description of the data, particularly in the region 
close to the maximal allowed rapidity. Some possibilities are discussed in the 
Appendix. 
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5 Summary and discussion 



Using the data on pseudorapidity distributions in D — Au collisions at 200 
GeV cm. energy j^j, we have shown that they can be reasonably well de- 
scribed by the wounded nucleon model [I]. This allows one to determine the 
contribution F(y) to particle production from one wounded nucleon which is 
a novel information, hitherto not available. The data show that (a) F(y) ex- 
tends over the full rapidity range, far beyond the hemisphere of the wounded 
nucleon in question, and (b) one observes a striking diference between the 
antisymmetric and symmetric parts of F(y). 

The last feature suggests that F(y) is built from two components, rep- 
resenting particle emission from two different sources. Extending the ideas 
formulated in [Hj, we proposed to identify these two sources as (i) the valence 
part of the nucleon and (ii) the soft part of the nucleon structure, dominated 
by gluons. This idea, accompanied with the assumption of the approximately 
flat gluon rapidity spectrum, explains immediately the striking linear behav- 
ior of the antisymmetric part of F(y), determined by the contribution from 
the valence source (the gluon contribution is symmetric in rapidity and thus 
does not contribute to the antisymmetric part of the spectrum). 

We thus conclude that the new data on D — Au collisions allowed us to 
obtain a qualitatively new information on particle production and to identify 
the two distinct sources inside the nucleon. 

Several comments are in order. 

(i) It should be emphasized that the model of wounded nucleons implies 
that the intensity of particle emission from a wounded nucleon does not 
depend on number of its interaction with the target. In our interpretation this 
means that the number of colour sources per unit of rapidity (in one nucleon) 
is independent of the number its interactions, i.e., independent of the number 
of colour exchanges between the projectile and target. The experimental 
verification of the model shows that such a saturation is indeed present. 

(ii) It seems likely that this saturation of particle emission is related to 
the concept of formation zone ^2] 5 i-e., strong reduction of soft emission 
from sources too close in rapidity. It would be interesting to investigate this 
question in more detail. 

(iii) Our interpretation of the data has very much in common with the 
dual parton model (DPM) in particular, our "valence" contribution 
corresponds to that of diquark-quark string in DPM, while our gluonic strings 
are analogous to the sea-quark strings of DPM. Ignoring the technical details 
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(inessential at this stage of discussion), the main difference is in the way we 
count the number of emitters. Although the number of valence sources is the 
same in the two models, counting of the "short" strings seems substantially 
different. In the dual parton model the number of the "short" strings equals 
the number of interactions between the projectile and target. As already 
explained above in (i), in our approach this is a property of the wounded 
nucleon, independent of the number of its interaction with the target. 

(iv) The simple distributions of partons and of emitted clusters used in 
our discussion were taken only for illustration. If more precise description 
of data is attempted, they must be accordingly modified, particularly in the 
region close to maximal rapidity. One example of possible modification is 
shown in the Appendix. It would be certainly interesting to perform such an 
analysis when the final version of data is available. 

(v) It was shown in that the wounded nucleon model does not describe 
correctly the data for the Au — Au collisions at RHIC energies. In particular, 
the particle density in the central rapidity region increases much faster than 
the number of wounded nucleons. It will be very interesting to compare these 
deviations with the ones observed in the present paper (c.f. Fig. 4). 



6 Appendix 



One possibility of the more adequate description of the data, still retaining 
the salient features of the model, is described below. Let 



p{y;y + ,y~ 



l _ e -(s/+-s/)A' 



h(y') 



e -(»-ir)A 

-(y-+Y)/n 



2Y 



v 



(21) 



One sees that these formulae modify the regions close to the maximal ra- 
pidities. In the limit of vanishing parameters /i, v and A one recovers the 
simple situation described in Section 4. When (??) is substituted into (??) 
one obtains 



dN(y) 
dy 



V>_ _ W(Y + y)W(Y -y) 
~ ~ 2Y - /x [1 - e~ 2Y /»] 



[22) 
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with 



W(z) =z-n 



and 



1 + -^—e'^ 
A — ji 



dN(y)W _ W(Y + y)V(Y-y) ]a .. 



-A 



1 + -*-e-'» 
ji — A 



Z7i (23) 



2F - /i [1 - e~ 2y /^] 



I6ii 



where 



V {z) = 1 + -^e- (Y -y^ + — l-e"^/^ 



v- A 



A-z/ 

For \ = jj, = u one obtains 

= ^ (l + e- 2/ ^) - 2/i (l - e-^) /8i 
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(24) 



(25) 



(26) 



V(z) = 1 - [1 + e- 2/ "/9i (27) 

These formulae show explicitly that the distributions are modified in a 
finite region close to the phase space boundary, where (\Y\ — \y\) is not too 
large. As one moves out from the boudary, i.e. when (\Y\ — \y\) are large, the 
corrections vanish exponentially, and we recover the results given by (??) and 
(??). The size of the region where the corrections are important is controlled 
by the parameters yU, A, and v. 

We have checked that these formulae are flexible enough to account for 
the results of Figs. 5 and 6. 
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